Quantum Hall conductance of two-terminal graphene devices 
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Measurement and theory of the two-terminal conductance of monolayer and bilayer graphene in 
the quantum Hall regime are compared. We examine features of conductance as a function of gate 
voltage that allow monolayer, bilayer, and gapped samples to be distinguished, including N-shaped 
distortions of quantum Hall plateaus and conductance peaks and dips at the charge neutrality 
point. Generally good agreement is found between measurement and theory. Possible origins of 
discrepancies are discussed. 



I. INTRODUCTION 

Graphene monolayers and bilayers are recently dis- 
covered two-dimensional gapless semimetals. The Dirac 
spectrum of excitations in monolayer graphene gives 
rise to a number of novel transport properties, includ- 
ing anomalous quantized Hall conductance with plateaus 
at 4(n +1/2) e^/fe, n = 0,±1,±2,... in multi-terminal 
samples .'321 Bilayer graphene has a quadratic, electron- 
hole-symmetric excitation spectrum, leading to quan- 
tized Hall conductance values Ane^/h, n = ±1,±2, ...W^ 
Both monolayer and bilayer graphene have a zeroth Lan- 
dau level, located at the charge neutrality point (CNP), 
which is eightfold degenerate in bilayers and fourfold de- 
generate in monolayers. Other Landau levels are all four- 
fold degenerate in both types of graphene.^ ^ The novel 
transport signatures not only reflect this underlying band 
structure, but serve as an experimental tool for iden- 
tifying the number of layers and characterizing sample 
quality.01 

In recent work on graphene, two-terminal magnetocon- 
ductance has e merged as one of the main tools of sample 
characterization.EEEl While a two-terminal measurement 
is not as straightforward to interpret as the corresponding 
multi-terminal measurement^, it is the simplest to per- 
form and may be the only measurement possible, for in- 
stance with very small samples. The presence of non-zero 
longitudinal conductivity causes quantum Hall plateaus 
measured in a two-terminal configuration to not be as 
well quantized as in multiprobe measurement.'^ As dis- 
cussed in detail below, plateaus exhibit a characteristic 
N-shaped distortion arising from the finite longitudinal 
conductivity that depends on device geometry. 

In this Article, we systematically examine two- 
terminal conductance in the QH regime for monolayer 
and bilayer graphene for a variety of sample aspect ra- 
tios (Table I). We especially focus on the features that 
can help to distinguish monolayer and bilayer graphene: 
the conductance extrema in the N-shaped distortions of 
the quantum Hall plateaus and at the CNP. We find that 
these features depend both on the sample aspect ratio 
and the number of graphene layers. Results are compared 



to recent theorjff^^ jn which two-terminal conductance for 
arbitrary shape is characterized by a single parameter ^, 
the effective device aspect ratio (^ = L/W for rectangu- 
lar samples, where L is the length or distance between 
contacts, and W is the device width). 

TABLE I: Measured two-terminal graphene devices. 



Sample Layers (Inferred) (L, W) [/xm] ^fit 



Al 


Monolayer 


(1.3, 1.8) 0.7 1.7 


A2 


Monolayer 


(0.4, 2.0) 0.2 0.2 


Bl 


Bilayer 


(2.5, 1.0) 2.5 0.8 


B2 


Bilayer 


(0.3, 1.8) 0.2 0.3 


C 


Monolayer 


Asymmetric 0.9" 0.9 



"effective aspect ratio, see Sect. IV. 

In Ref. |lll . the positions of conductance extrema on 
the distorted plateaus were linked to incompressible den- 
sities. Here we find that this relation can be used to 
distinguish monolayer and bilayer graphene devices even 
when the distortions of the plateaus are strong. The anal- 
ysis of rectangular two-terminal samples is extended to a 
sample with asymmetric contacts, extracting an effective 
sample aspect ratio via conformal mapping. Best-fit val- 
ues of the aspect ratio, ^gt , obtained by fitting the theory 
to the experimental data, are compared to the measured 
sample aspect ratio, ^g. Agreement is generally good, 
but not uniformly so. We speculate on possible causes 
of these discrepancies, including inhomogeneous contact 
resistance, electron and hole puddles, and contributions 
of transport along p-n interfaces. 

A. Qualitative Discussion 

Representative theoretical plots of two-terminal con- 
ductance for monolayer, bilayer, and gapped bilayer 
graphene as a function of filling factor, v, are shown in 
Fig. 1. For both monolayers and bilayers, the absence of 
an energy gap between the conduction and valence bands 
gives rise to a zero energy Landau level (LL^, which can 



2 



either increase or decrease the two-terminal conductance 
around the charge neutrality point, depending on the as- 
pect ratio of the sample. The eightfold degeneracy of the 
zero-energy LL in bilayer graphene^ enhances the size of 
this feature relative to monolayer graphene. 

A gap in the spectrum of bilayer graphene opens when 
the on-site energy in one layer differs from the on-site 
energy in the other P^l This may result, for instance, from 
asymmetric chemical dopin^Sl or electrostatic gating)^ 
The gap splits the zero-energy LL, suppressing conduc- 
tance at the CNP. The qualitative effect of a gap in the 
bilayer spectrum can be seen in Fig. 1 by comparing the 
gapped case [Fig. 1(c)], which always has a zero of con- 
ductance a,t — 0, to the gapless cases [Figs. l(a,b)], 
which has a non-zero value of conductance a,t v = 0. 

Also illustrated in Fig. 1 is how the aspect ratio of the 
sample affects the two-terminal conductance near quan- 
tum Hall plateaus for all three spectrum types. Finite 
longitudinal conductivity leads to N-shaped distortions 
of the plateaus,'Sl which are of opposite signs for aspect 
ratios ^ < 1 and ^ > 1. Note, however, that the ex- 
trema of conductance — minima for ^ < 1 and maxima 
for ^ > 1 — are aligned with the plateaus centers, which 
coincide with the incompressible density values (different 
for monolayers and bilayers). Distorted plateaus thus re- 
main useful for characterizing the number of layers and 
density. 

The back-gate dependence of conductance for the five 
samples reported are most similar to those in Figs. l(a,b), 
indicating that these samples are single layers and gap- 
less bilayers only (see Table 1). We use the model of 
Ref. [IT] to fit the conductance data treating the aspect 
ratio as a fit parameter. In doing so, our presumption is 
that the visible dimensions of the sample may not reflect 
the actual pattern of current flow. Since the conductance 
problem for a sample of any shape can be reduce d to that 
of an effective rectangle via a conformal mappin^^^^^^, 
which depends on the sample shape but not on the con- 
ductivity tensor, the rectangular geometry is universal for 
two-terminal conductance. Thus the model of a conduct- 
ing rectangle with an unspecified aspect ratio is suitable 
for describing systems in which current pattern is not 
precisely known. 



B. Sample Fabrication and Measurement 

Graphene devices were fabricated by mechanically ex- 
foliating highly oriented pyrolytic graphite^'^ onto a n"'"'" 
Si wafer capped with ~ 300 nm of Si02. Potential single 
and bilayer graphene flakes were identified by optical mi- 
croscopy. Source and drain contacts, defined by electron 
beam lithography, were deposited by thermally evapo- 
rating 5/40 nm of Ti/Au. The aspect ratio, ^s, of each 
sample was measured using either optical or scanning 
electron microscopy. 

Devices were measured in a ^He refrigerator allowing 
dc transport measurements in a magnetic field |_B| < 8T 
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FIG. 1: (color online) Theoretical two-terminal QH conductance 
g as a function of filling factor u (Ref. ,11| ) shown for (a) single- 
layer graphene (b) bilayer graphene (c) gapped bilayer graphene 
for effective aspect ratios § = L/W = 2 (black) and 0.5 (red). 
Finite longitudinal conductivity due to the states in the middle of 
each Landau level distorts the plateaus into N-shaped structures, 
which are of opposite sign for 5 < 1 and ^ > 1. Local extrema 
of g at filling factors u = ±2, ±6, ±10... for single layers and at 
u = ±4, ±8, ±12... for bilayers are either all maxima (§ < 1) or all 
minima (^ > 1). For gapless monolayer and bilayer samples (a,b), 
glu = 0) is a minimum for ^ < 1 and maximum for ^ > 1; for the 
gapped bilayer (c) g vanishes at = for all 5. 

perpendicular to the graphene plane. Unless otherwise 
noted, all measurements were taken at base temperature, 
T ~ 250mK. Differential conductance g = dl/dV, where 
/ is the current and V the source-drain voltage, was mea- 
sured using a current bias (/ chosen to keep eV < ksT) 
and standard lock-in technique at a frequency of 93Hz. 
All samples show B — Q characteristics of high-quality 
single-layer and bilayer graphen^^^J a CNP positioned 
at back-gate voltage Vbg ^ and a large change in g (in 
excess of 20e^//i) over the Vbg range of ±40 V. 

II. MONOLAYER SAMPLES 

Figure 2(a) shows the two-terminal conductance ^(Vbg) 
for sample Al (^^ = 0.7) at B= 8 T (black trace). 
Plateaus are seen at = ±2 near — but not equal to — 
2e^/h, with values of ~ 2.3(2.7) e^//i on the electron 
(hole) side of the CNP. At the CNP (Vbg - 2.3V), g 
departs from the quantized values, dropping to a mini- 
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mum of ~ lAe^/h. At higher densities, the conductance 
exhibits a series of maxima with values sHghtly above 6, 
10, lAe^/h. Maxima on the hole side consistently have 
slightly higher values, a feature observed in all the sam- 
ples measured. The inset of Fig. 2(a) shows g in the 
QH regime as a function of Vbg and B. Dashed black 
lines indicate the filling factors = Ugh/ eB (where Ug is 
the carrier density) of —6, —10, and —14 and lines align 
with the local maxima of g{Yog,B). Vbg was converted 
to Ug using a parallel plate capacitance modeP^, giving 
Ug = Q!(Vbg + Kffsot) with a = 6.7 X lO^^cm-^V'^ and 

^offset = 2V. 

Measured ^(Vbg) [black curve in Fig. 2(b)] for sam- 
ple A2 (^s = 0.2), made using the same graphene flake 
as Al, shows distinctive differences from the measured 
(7(Vbg) of sample Al. In particular, at the CNP (l^g — 
— 1.5V), g exhibits a sharp peak with a maximal value 
of ~ 8.8e'^/h. Away from the CNP, the conductance has 
maxima which are much stronger than those of sample 
Al. The inset of Fig. 2(b) shows g(Vbg,S). For this 
sample, the dashed lines representing the incompressible 
filling factors ±6, ±10, ±14 now align with the minima 
in g. Here we used the Vbg to ng conversion factors of 
a = 6.7 X 10"'^°cm~^V~^ (the same as for sample Al) and 

Offset = -I.IV. 

The observed features in g for samples Al and A2 can 
be compared to theorjOT for two-terminal quantum Hall 
conductance, which uses a model of a conducting rect- 
angle L X W with a spatially uniform conductivity. The 
filling factor dependence of the conductivity tensor is ob- 
tained using the semicircle relation for quantum Hall sys- 
tems, derived in Ref. [19', which is applied independently 
for each Landau level. Landau level broadening due to 
disorder is included in the theory as a gaussian broaden- 
ing e~-^^'^~'^"'> , where Vn is the center of the LL and A 
is a fitting parameter. The total conductivity tensor is 
taken to be a sum of the contributions of individual Lan- 
dau levels. The current-density distribution for a rect- 
angular sample with an arbitrary aspect ratio is found 
analytically by conformal mapping (see Refs. |15I16I17] ). 
The current density is then integrated numerically along 
suitably chosen contours to evaluate total current and 
voltage drop, from which g — I/V is obtained. 

Along with the experimental traces, Figs. 2(a,b) also 
show the theoretical curves for the best-fit (solid red 
trace) and the actual sample aspect (dashed blue trace) 
ratios. For sample Al, the best-fit aspect ratio, ^gt = 1.7, 
differs considerably from the measured value, S,g — 0.7. 
For sample Al, the best fit gives the Landau level broad- 
ening parameter A = 1.2. This theoretical curve (ffit — 
1.7) reproduces the essential features of the data: local 
maxima align with the filling factors ±2, ±6, ±10, and 
g exhibits a dip at the CNP. 

The alignment of conductance minima with densities 
corresponding to the integer filling factors as well as a 
peak at the CNP observed for sample A2 are consistent 
with theoretical predictions for a short, wide monolayer 
graphene sample. As illustrated in Fig. 2(b), the best-fit 
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FIG. 2: (color online) (a) Inset: Conductance g in the quan- 
tum Hall regime as a function of B and Vbg at T = 250 mK for 
sample Al. Black dashed lines correspond to filling factors of 
u = —6, —10, —14, —18 and align with the local maxima of con- 
ductance. Main: (black) Horizontal cut of inset giving g{V\,g) at 
B = 8T and calculated g for the best-fit equivalent aspect ra- 
tio ^Ht = 1-7 (solid red curve) and the actual sample aspect ratio 
5s = 0.7 (dashed blue curve) using Landau level broadening param- 
eter A = 1.2. (b) Inset: Conductance g in the quantum Hall regime 
as a function of B and Vbg at T = 250 mK for sample A2. Black 
dashed lines correspond to filling factors of f = —6, —10, —14, —18 
and align with the local minima of conductance. Main: (black) 
Horizontal cut of inset giving g(Vbg) at i? = 8T and calculated g 
ior 5flt = 0.2 (solid red curve) and = 0.2 (dashed blue curve) 
(A = 1.2, the same as sample Al). 



aspect ratio ^gt = 0.2 agrees well with the measured £,g 
for sample A2. 

We observe that the size of peaks and dips in 
Fig. 2(a,b) increases for higher LL. In contrast, theorjE] 
predicts that peaks and dips at > LLs are all roughly 
the same. This discrepancy may reflect the inapplicabil- 
ity of the two-phase model approach of Ref. [12] , which 
underlies the semicircle law obtained in this work, to 
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FIG. 3: (color online) Inset: Measured g of sample Bl as a func- 
tion of B and Vbg at T = 4K. Black dashed lines, corresponding 
to u = —12,— 16, —20, align with local minima of g. No minima 
are observed at i/ = 8 for 5T < B < 8T. Main: Horizontal cut of 
inset at _B = 8T (black), and calculated g using A = 0.7 for = 
2.5 (dashed blue curve) and = 0.8 (solid red curve). 



FIG. 4: (color online) Measured g(Vbg) for sample B2 (black) and 
the calculated g using A = 0.25 for = 0.2 (dashed blue trace) and 
^fit =0.3 (solid red trace). Two key features in the curve suggest 
this sample is a gapless bilayer, namely, a pronounced peak in g 
near the CNP, and the larger spacing between the two minima 
straddling the CNP compared to the spacing (5Vbg ~ 9.5 between 
other consecutive minima. 



higher LLs. Indeed, because for Dirac particles the spac- 
ing between LLs decreases at higher energies as an in- 
verse square root of the level number, one may expect 
mixing between non-nearest LLs to increase at high en- 
ergies. Such mixing can lead to the longitudinal conduc- 
tivity values in excess of those of Ref. [19] , which only 
considers mixing between nearest LLs (see the discussion 
in Ref. [ID]). 

To take these effects into account, we extend the model 
of Ref. [11] by assuming that the contribution of the n^^ 
LL to the conductivity tensor in monolayer graphene is 
described by a modified semicircle ("elliptic") law, 

SnCrl^ + Al{Sn<7xy - (^ly ,n){^nCr xy - Cr°y^„,) = 0, (1) 

where 5n<Jxx and Snf^xy are the effective longitudinal and 
Hall conductivities, cr^y „ and are the quantized 

Hall conductivities at the neighboring plateaus. Here n 
and n' are neighboring LL indices, related by n' = n + 1 
(except the doubly degenerate v =0 LL for the bilayer, 
in which case n — —1 and n' — 1). The An account for 
departures from the semicircle law. We take An ~ 1 for 
n = 0,±1, and An ~ 2 for other LLs, consistent with 
previous observations.!^ 

III. BILAYER SAMPLES 

The black curve in Fig. 3 shows measured g{Vhg) for 
sample Bl (^^ = 2.5) at B = 8T and T = 4K. This 
sample has two features indicating that it is a bilayer 
sample: plateaus in conductance appearing near 4, 8, 12 



and 16e^//i, and a conductance maximum at the CNP 
whose relative size is much larger then those at higher 
LLs. The conductance values at the plateaus v = ±4 here 
are lower than the expected Ae^/h for a bilayer sample, 
falling to 2.7(3.1) /h on the electron (hole) side of the 
CNP. The peak value in conductance at = (Vbg — 
0.5 V) is 5e^/h. At higher filling factors, the plateaus 
exhibit two different behaviors, showing a flat plateau at 
1/ = 8 and a plateau followed by a dip at h' — 12. The 
small dips align with the filling factors v — —12, —16, —20 
for 5T < B < ST (see inset of Fig. 3), using a = 7.2 x 
10i°cm-2v-i and VoSsct = 0.5V. 

Theoretical g curves for aspect ratios £,s = 2.5 (dashed 
blue curve) and = 0.8 (solid red curve) are shown 
in Fig. 3. Theoretical 5 (Vbg) curves for these two aspect 
ratios are similar at high density, but differ for v — 0: the 
curve for = 2.5 has a dip in conductance at the CNP 
while the best-fit curve (^gt = 0.8) has a peak, similar 
to the experimental curve. The curve for = 0.8 also 
agrees better with experiment at higher densities. 

In some cases the quantized conductance values are 
found to be quite different from the expected quan- 
tized values, as demonstrated in Fig. 4 for sample B2 
(^s = 0.2). In this sample, g reaches a maximum of 
13.5 eV^ at the CNP, with adjacent minima of 5e^/h. 
Away from the CNP, conductance plateaus appear at val- 
ues of ~ 16e^/h and 23e^/h, neither of which are near 
expected values for monolayer or bilayer graphene. Since 
there are no strong peaks or dips in g away from charge 
neutrality, as is expected for a device with a <C 1, it is 
difficult to determine the number of layers from the loca- 



5 




FIG. 5: (color online) Measured g{Vhg) for sample C (black) 
and calculated conductance (solid red curve) for = 0.9 (A = 
0.7). The asymmetric contacts of this sample can be conformally 
mapped onto a rectangle, producing a device aspect ratio of t;s = 
0.9 (dashed red curve, directly under solid red curve). 



tion of the conductance extrema. There are two conduc- 
tance features, however, that suggest the sample is gap- 
less bilayer graphene. First, the peak at = is much 
more pronounced than any other peak in the conduc- 
tance. Second, the spacing in Vbg between the two lowest 
LLs is twice as large as the spacing between any other two 
successive LLs (in Fig. 4, SVhg ~ 9.5 V). Both features 
arise in bilayers as a result of the zero-energy LL being 
eightfold degenerate, twice as much as all other bilayer 
LLs and the zero-energy LL in single layer graphene.^ 
The theoretical prediction, using = 0.3 (solid red line) 
and = 0.2 (dashed blue line), for the bilayer sample 
B2 are shown in Fig. 4. 



IV. NON-RECTANGULAR SAMPLES 

In this section we extend the comparison of theory and 
experiment to a non-rectangular device, sample C, shown 
schematically in the inset of Fig. 5. The measured two- 
terminal conductance of sample C (black curve in Fig. 5) 
has properties very similar to those expected for a square 
monolayer sample: around the CNF the conductance is 
nearly flat with value close to 2e^//i, monotonically in- 
creasing on the electron and hole sides at filling factors 

Theoretical curve shown in Fig. 5 is obtained from the 
conducting rectangle model using a best-fit effective as- 
pect ratio ^fit = 0.9 and the LL broadening parame- 
ter A = 0.7. This choice of parameters yields particu- 
larly good agreement for < 6. At higher fillings, the 
plateaus are washed out, suggesting that the LL broad- 
ening is stronger for LLs |n| > 2. It is interesting to 
compare the best-fit value to an effective aspect ra- 
tio, obtained from conformal mapping of sample C to 
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FIG. 6: (color online) A polygon representing sample C (see 
Fig. 5). Black regions correspond to contacts, length scale 
a = 200 nm. 

a rectangle. As discussed below, this conformal map- 
ping can be constructed directly, owing to the relatively 
simple geometry of sample C. The effective aspect ratio 
obtained in this way is ~ 0.9, which is consistent with 
the best-fit value. 

Before we proceed to construct the conformal map- 
ping we note that the geometry of sample C, pictured in 
Fig. 6, is that of a polygon. In principle, any polygon 
can be mapped onto the upper half-plane by inverting a 
Schwarz-Christoffel mapping. However, since this map- 
ping is defined by a contour integral, the inverse mapping 
can only be found numerically. In order to circumvent 
this difficulty, two approximations are employed below, 
allowing the desired mapping to be constructed as a com- 
position of a few simple mappings. 

The steps involved in this construction are illustrated 
in Fig. 7. First, the rectangular shape in Fig. 6 is re- 
placed by a semi-infinite strip shown in Fig. 7(a). This 
approximation should not significantly affect the conduc- 
tance, as the current flows mostly in the region between 
contacts 1-2 and 3-4. Without loss of generality we set 
the length scale a — 1. 

Our next step is to straighten out the contact 3-5-6-4. 
For that, let us consider an auxiliary mapping that maps 
the upper w plane onto the upper z plane with a removed 
rectang 

We choose the parameter A to be equal 

so that the removed rectangle has vertices 

^3^4 = ±A, Z5,6 ^±A + iA. (4) 

These points correspond to the points w^^i — ±\/2, 
ws^e — il in the w plane. The value of A ensures that 
the edge of the sample on the x axis remains on the x 
axis under the mapping ([2| . The distance between points 
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FIG. 7: (color online) Three steps used to map the polygon 
in Fig. 6 (sample C) onto the upper half-plane (schematic). 
First, the rectangle in Fig. 6 is replaced by a half-infinite 
strip, extending indefinitely to the right (a). Next, we map 
the domain shown in (a) onto a rectangle with contact 3-5- 
6-4 straightened out (b). Under this mapping, the sample is 
slightly distorted, as indicated by the grey polygon in (b). 
Because the deviation of the grey polygon boundary from the 
original sample boundary (red line in (b)) is fairly small, it 
can be neglected, giving a half-infinite strip (c). Finally, the 
domain (c) is mapped onto the upper half-plane (d), which 
allows to find the cross ratio A1234, Eq. ([9|, and evaluate the 



effective aspect ratio, Eq. ( 10 1 



^3 and ^5 plane equals A, as follows from Eq. (|2| and the 
identity 



75 



^2 _2 



1/2 



e~2 



1/2 



d^, (5) 



which can be proved by making the change of variables, 
= y/2 — X in the integral in the left-hand side of Eq. ^ , 
and ^ — y/x in the integral in the right-hand side of 
Eq. (|5|. 

The removed rectangle has aspect ratio equal to 2, the 
same as that for the contact 3-5-6-4, however, their di- 
mensions differ by a factor of A. Scaling and shifting 
both z in w. 



z — A{z — 5), w — A{w — 5), 



(6) 



we obtain the required mapping which straightens out 
the contact 3-5-6-4. 

The second approximation is necessary because the 
mapping ([2| , ([6]) , while straightening the segments 3-5-6- 
4, distorts the rest of the boundary. We notice, however, 
that sufficiently far from the contact 3-5-6-4 the mapping 
([2]) is close to the identity: 



z(w > 1) = w -I- 0(l/w), |z-5|>l. 



(7) 



This property and the relatively small size of the seg- 
ments 3-5-6-4 compared to the strip width guarantees 
that the distortion is small. This is shown schematically 
in Fig. 7(b), where the curved grey polygon represents 



the actual image of the sample, with the deviation of its 
boundary from the strip of the same asymptotic width 
(shown in red and exaggerated for clarity). The devia- 
tion is indeed small: by investigating the mapping ^ , ^ 
numerically we found that the boundary is displaced the 
most at the point 2 which is shifted by approximately 0.3 
away from its original position 2' along the real axis. This 
is small compared to the sample width, equal to 6, which 
allows us to neglect the displacement of the boundary. 
Thus we assume that the mapping ([2|, ^ transforms 
sample C into the semi-infinite strip shown in Fig. 7(c). 

After this approximation is made, it is straightforward 
to transform the semi-infinite strip in Fig. 7(c) into the 
upper half-plane, which can be done by the following 
mapping. 



C — cosh- 



6 



(8) 



In the C plane, the contacts are mapped on the real axis, 
with the end points 1, 2, 3 and 4 mapped to Ci = C2 = 
1, C3 ~ 2.11, C4 « 23.57. From these values, following the 
procedure described in Ref . |llj (Appendix) , we compute 
the cross ratio 



il234 



(Cl-C4)(C3-C2) 
(Cl-C2)(C3-C4) 



-0.64, 



(9) 



and then obtain the aspect ratio from the relations 



6- 



L _ K{k') 
W ~ 2K{ky 



A1234 - (1 - fc')/2fc, 



(10) 



where K{k) is the complete elliptic integral of the first 
kind, and fc' = (1 — A:^)^/^. This procedure yields the 
value = 0.9, identical to that found from the best fit 
to a conducting rectangle model (see Fig. 5). 



V. SUMMARY AND DISCUSSION 

In summary, we have studied the effect of geometry 
on the conductance of two-terminal graphene devices in 
the QH regime, comparing experiment and theory. The 
quantized QH plateaus typically exhibit conductance ex- 
trema that are stronger for wide, short samples. For non- 
rectangular samples, the equivalent rectangle approach 
appears works well. 

Theoretically, for short, wide samples (^^t < 1) the 
two-terminal conductance a conductance minimum is ex- 
pected at filling factors where plateaus exist in multiter- 
minal devices, while for long, narrow samples (^gt > 1), 
a conductance maximum is expected at these filling fac- 
tors. Along with the behavior at the CNF, these sig- 
natures provide a clear way to identify the number of 
layers in the sample even when the quantization is weak 
or absent. 

We find in the five samples measured that conductance 
as a function of gate voltage is well described by theory, 
allowing for a phenomenological Landau Level broaden- 
ing, and treating the aspect ratio as a fit parameter. In 
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some samples, however, the best fit aspect ratio differs 
considerably from the measured aspect ratio of the sam- 
ple. 

What might lead to the discrepancy between some of 
the measured and fit aspect ratios? One source of dis- 
crepancy could be that only part of the contact actu- 
ally injects current. It might also be that the contacts 
locally dope the graphene, causing the aspect ratio to 
appear smaller. This latter scenario, however, would 
require that the doping penetrates ~ 500 nm into the 
graphene, which is ~ 2 orders of magnitude larger than 
expected.l^^l Another, more interesting possibility could 
be that the picture of an effective medium characterized 
by local conduction, on which the argument leading up 
to the semi-circle relation^ is based, may not hold. This 
might arise, for instance, from large density fluctuations, 



giving rise to electron and hole puddlei^ forming a net- 
work of p-n interfaces along which conduction occurs. In 
this case, the effect of the back gate is to alter the per- 
colation properties of this p-n network. Transport me- 
diated by such states would almost certainly change the 
conventional picture of local conduction. Further studies 
are required to clarify the physical mechanism responsi- 
ble for the observed behavior. 
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